Long quantum channels for high-quality entanglement transfer 
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High-quality quantum-state and entanglement transfer can be achieved in an unmodulated spin 
bus operating in the ballistic regime, which occurs when the endpoint qubits A and B are coupled 
to the chain by an exchange interaction jo comparable with the intrachain exchange. Indeed, the 
transition amplitude characterizing the transfer quality exhibits a maximum for a finite optimal 
value jQ^*(A^), where N is the channel length. We show that jQ^^(A^) scales as A^~-^/^ for large 
N and that it ensures a high-quality entanglement transfer even in the limit of arbitrarily long 
channels, almost independently of the channel initialization. For instance, the average quantum- 
state transmission fidelity exceeds 90 % for any chain length. We emphasize that, taking the reverse 
point of view, should jo be experimentally constrained, high-quality transfer can still be obtained 
by adjusting the channel length to its optimal value. 
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I. INTRODUCTION 

Quantum state transfer between distant qubits (say, 
A and B) is a fundamental tool for processing quantum 
information. The task of covering relatively large dis- 
tances between the elements of a quantum computer, 
much larger than the qubit interaction range, can be 
achieved by means of a suitable communication channel 
connecting the qubits A and B. 

Spin chains are among the most studied channel pro- 
totypes [1]. In particular, the XY S = l/2 model has 
been widely employed as a tool for testing and analyzing 
quantum communication protocols. For such a commu- 
nication channel to be experimentally feasible one could 
be lead to consider systems with uniform intrachannel 
interactions [2-8] and whose operation does not require 
peculiar initialization procedures, although other pro- 
posals have been recently put forward |9l-[l2]. On the 
other hand, the quantum-transfer capabilities of homo- 
geneous spin channels have not yet been fully explored 
under many respects. For instance, it is a common belief 
that the longer the chain the worse is the transmission 
fidelity [ll,[2|,ll3l,ITl as an effect of dispersion, and chains 
up to only a few tenths of spins are often considered. 

The present work is devoted to the study of the ballistic 
regime, where the transmission can be depicted [6, 15, 16] 
in terms of a traveling wavepacket carrying the infor- 
mation about the state of the endpoint qubit A, even- 
tually yielding the state reconstruction at the opposite 
endpoint qubit B thanks to the overall system's mir- 
ror symmetry [ol-llll, fl7| . The ballistic regime differs 
from that arising in the limit of very weak endpoint cou- 
plings d, m, [H, [iQ^ where (almost) perfect state transfer 
occurs at very long times as a result of a Rabi-like pop- 
ulation transfer involving only the lowest-energy single- 
particle modes. Understanding the basic mechanism of 
ballistic transfer, where the number of involved single- 
particle modes will be shown to be of the order of A/"^/^, 
allows us to devise an optimal value of the endpoint in- 



teractions for any A", and vice versa. Remarkably, the 
corresponding transmission quality, as witnessed by the 
state- and the entanglement fidelity, does not decrease to 
zero when the channel becomes very long, but remains 
surprisingly high. 

We consider the setup illustrated in Fig.[TJ the channel 
connecting the qubits A and B is a one-dimensional array 
of N localized S = 1/2 spins with exchange interactions 
of XX Heisenberg type and a possible external magnetic 
field applied along the z direction. This gives the total 
Hamiltonian the following structure 
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Y,{sfsf^,^sysy^^-h^s^ (i) 
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-joE {SfSUi+SfS!^,)-ho{S§+S^^+,) , 

i=0,N 

where the qubits A and B sit at the endpoint sites 
and A'+l of a one-dimensional discrete lattice on whose 
sites 1, 2, A" the spin chain is set. The exchange inter- 
action (chosen as energy unit) and the magnetic field h 
are homogeneous along the chain, and an overall mirror 
symmetry is assumed, implying the endpoint coupling jo 
and field ho to be the same for both ends. The N spins 
constituting the XX channel are collectively indicated 
by r. We will focus our attention on how the state of the 
qubit B evolves under the influence of the chain F, and 
depending on the initial state of the qubit A; the latter 
is possibly entangled with an ancillary qubit C. The re- 
sults of the analysis are used for gathering insights on the 
quantum-information transmission through the chain, so 
as to characterize the dynamical evolution of the overall 
system and to maximize the quality of the quantum-state 
transfer. 

Even though the overall scheme could also be used for 
realizing tasks other than quantum information transfer, 
via the dynamical correlations that the chain induces be- 
tween A and B [20], our approach is specifically tailored 
for studying transfer processes along the chain: the qubit 
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FIG. 1. The endpoints of a quantum channel V are coupled 
to the qubits A and B, via a tunable interaction jo; A can be 
entangled with an external qubit C. 



density operator this can always be expressed as 

a linear function of a suitable input density matrix p^", 
which may be the initial one of qubit B itself or that of 
any other qubit playing a role in its evolution: 

p(t) = St r , (2) 

where the quantum operation St is a trace-preserving, 
completely positive, convex-linear map from input den- 
sity operators p^^ to density operators of the qubit at 
time t (see, e.g., Ref. Hil). The effect of the linear map 
can be represented in terms of a 4 x 4 time-dependent 
matrix 



B or the qubit-pair BC are considered as target system, 
depending on whether the quantum-state of the qubit A 
or that of the qubit-pair AC are to be transferred, re- 
spectively. 

In Section [iT] the formahsm used to study the dynami- 
cal evolution of the composite system is introduced and 
we derive the corresponding time-dependent expressions 
for the quantities used for estimating the quality of the 
quantum-state and entanglement transfer processes. In 
Section inil the proposed formahsm is applied to the XX 
model described by Eq. ([T]). In Section |IV] we put forward 
an analytical framework in order to improve the under- 
standing of the conditions inducing an optimal ballistic 
dynamics. The resulting high-quality transfer processes 
along the spin chain are analyzed in Section |Vl Conclu- 
sions are drawn in Section IVH where comments about 
possible implementation of the procedure are also put 
forward. Relevant details of calculations are reported in 
the Appendices. 



II. DYNAMICS AND TRANSFER QUALITY 
A. Dynamics 

The density matrix of a qubit B can be written as 
p = C/i where p = 1,...,4, are complex num- 

bers, and {C/i} is an orthonormal basis in the Hilbert 
space of 2 X 2 matrices endowed with the Hilbert-Schmidt 
trace-product, Tr[C^Ci^] = (5^^^. In some cases it is useful 
to choose Hermitian C^, for instance the identity and the 
three Pauli matrices. However, a more direct connection 
with the computational basis {|0) = It)^ |1) = foi" 
the pure states of B is established by choosing = 
and setting p = l-|-2iH-j, with z, j = 0, 1. In what fol- 
lows we will use this description, and hence use indices 
/i, z^, A running from 1 to 4, understanding summation 
over any repeated index, e.g., p = b^(^. Hermitic- 
ity of p implies hi and 64 to be real and 62 = 63, while 
Trp=l means 61 + 64 = 1, so that only three real pa- 
rameters are independent. The positivity of p requires 
that I62P <6i64. 

The dynamics of B is described by its time-dependent 



(3) 



In fact, by writing p{t) = h^{t) C,^ and /O™ = a^C^) one 
gets 



(4) 



If B is part of a larger system S ruled by a total Hamil- 
tonian H and prepared in the initial state it is also 



(5) 



where ( • ) = Tr[ • and (^{t) = e'^' IgxB ^ Cm e 

acting on the Hilbert space of the qubit B. The matrix 
elements T^u{t) explicitly follow from Eq. (j4]) and ([5j). 

Here we will essentially focus on the time evolution 
of the qubit B and take the input density operator in 
Eq. ([2]) to represent the initial state of the qubit A, as 
this is the most natural setup for studying quantum state 
and entanglement transmission processes from A to B. 
Referring to the specific setup described in Fig. [1] we 
prepare the overall system AUFUB in the initial state 



p^ p" Q<) p^ and let it evolve into e ^^^ptot^iHt^ 



where 1-L is the total Hamiltonian ([T]). Notice that 
has a fully separable structure: this is not a necessary 
condition for the explicit determination of the dynamical 
matrix (a structure such as p^ (g) p^^ would also work), 
but it is quite meaningful in our scheme, as the qubit B 
is assumed not to interact with the chain during the ini- 
tialization process, which makes quite artificial its being 
entangled with the chain before the dynamics starts. 

The density matrix of B evolves according to p^ (t) = 
TrAur [e-'^^p^^'^e'^^], from which 



T^At) = Tib [c^ TrAur [e"*^'C. ' 



(6) 



where ( • )iy = Tt[- ^ p^ ^ p^]. Notice that the initial 
state of the qubit A does not enter the expression of 
T^iy{t), setting this procedure in the framework of general 
tomographic approaches. 

Let us now consider the time evolution of the quan- 
tum state describing the qubit-pair CUB when the to- 
tal system is initially prepared in the state = 
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p^^ p^. By a trivial generalization of the descrip- 
tion introduced at the beginning of this section, we write 
P^^ =9^uC^^G^ and find 



This equation shows that T^u{t) is the only ingredient 
needed not only for deriving the dynamics of the qubit 
B, but also of the qubit pair CUB, providing the pair 
CUA was prepared separately from the rest of the system, 
and C is non- interacting [22]. In particular, if C and 
A are initially prepared in one of the Bell states, say 
(1 00) + |11))/V^, it is g^^ = ^(5^^^ (similar conditions 
hold for the other Bell states), i.e.. 



PBell 



Cm 



and hence 



(8) 



(9) 



The trace-preserving property of the dynamical map, 
Tr [ft C/i] = Tr together with the specific choice of 
the basis {C/i}, entails T12 + T42 = T13 + T43 = and 
Til + T41 = Ti4 + T44 = 1; the Hermiticity of p^ further 
implies that T^ia^ and T^/^a^ are real and T^2<^/i = T^^a"^^. 
Moreover, from Eq. ([6j) it clearly follows that symmetry 
properties of the total Hamiltonian and of the initial state 
p^ p^ may result in constraints for the matrix elements 
In particular, when both are symmetric under ro- 



tations around the z-axis, it is 



^^2 



T2ll = T22 S, 
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f /i3 



/i3 



(10) 



and only three matrix elements, say Tn, T22, and T44, 
need to be determined. 



B. Quality of transfer processes 

In order to study the quality of the transfer processes 
mediated by the spin chain, we specifically consider the 
entanglement transfer from CuA to CUB, when the for- 
mer spin pair is initially prepared in a pure state. In this 
case, the entanglement fidelity, which is a proper mea- 
sure of the quality of entanglement transmission, reads 
Tent{t) = Tr [pCAV^^(^)], which is, via Eq. 



(11) 



The entanglement fidelity measures how close the state of 
CUB at time t is to the initial state of CUA. In particular, 
if C and A are initially prepared in a Bell state, Eq. (j8j), 
it is 



x-Bell/ 



entit) = iT^^it) , (12) 

while the entanglement of the pair CUB is measured by 
C?|"W = |C[T.^(t)Cp®C.] , (13) 



where C[p] is the concurrence [23[ between two qubits 
in the state p. Note that Eqs. ©, ([I2|), and are 
independent of which Bell state is chosen, since different 
Bell states are connected to (|00) + |11)) /V2 by unitary 
operations on C, which is isolated. 

Another relevant tool for evaluating the quality of the 
transfer processes is the fidelity of transmission from A 
to B, which reads ^ab(^) = Tr [p^^ p^{t)] provided that 
p^ is a pure state, i.e., by Eq. (j4]), 

J',At) = b^{t)a; = T^,{t)a;a, . (14) 
If A is initially prepared in the state |V^6>(^) = cos f |0) + 



e^'^sinlll), meaning ai 



sin f , and a2 



Eq. (p!4|) can be averaged over all 
possible initial pure states by integrating over the Bloch 
sphere, resulting in 



6 : 



(15) 



which can be compared with Eq. (fT2]) to obtain the rela- 
tion H m 



3 -^ent ' 



(16) 



It is worth noticing that a high average fidelity could 
still allow for states that are poorly (or even not at all) 
transferred, while the ultimate goal is the transmission 
of any state. 

When the setup is such that Eqs. (p!Q]) hold, from 
Eqs. ([ni) and ([l3|) it follows that 

^fnfW = i(^ll+^44 + 2ReT22) , (17) 

C^f{t) = max{0, IT22I - . (18) 

Moreover, if A is prepared in the pure state \^l^e^)i 
Eq. (dH) explicitly reads 

-^AB W = ^11 cos^ f + ^44 sin^ f 

+ [2(1 + ReT22) - (Til + T44)] sin^ f cos^ f , (19) 

which shows that fo=Tii and /i=T44 are the trans- 
mission fidelities of the states |0) = IV^o^^) and |1) = |'07r(^), 
respectively, while / = ^(1+ ReT22) represents the trans- 
mission fidelity of the states (|0) + e*^|l))/V^= {ip^^). 
The fidelity (fTQ]) only depends on ^ G [0, tt] and its ex- 
trema can be easily determined. To this purpose rewrite 
Eq. (dl 



^AB W = / + Ufo-fi) COS + i(/o+/i - 2/) cos^ ; 

(20) 

it follows that ^ab(^) takes the values /o and fi at the 
extrema of the range cos^ G [—1,1] and can have a mini- 
mum 



/m — / 



1 ifo-fi? 
8 /0+/1-2/ 



(21) 
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/i-/o 



in between provided that Cm — 2(/^+/o-2/) 
Therefore the bounds of J^^^ (t) are: 



G (-1, 1). Eqs. (do]) hold, we find 



_ J min{/i,/o} if |ci„| > 1 

^'~\fm if |C„,|<1 

J-™f^(t)=max{/i,/o} . 



(22) 
(23) 



One can see that if /i = /o the best transmitted states are 
|0) and |1) which, on the other hand, are the best and 
the worst transmitted ones (or viceversa) if |cm| > 1. This 
pecuhar role of the computational states is obviously a 
consequence of the assumed symmetry. 



III. THE XX MODEL 
A. Dynamical evolution 

In this section we specifically consider the Hamilto- 
nian ([T]). The system AUFUB is prepared in the state 
ptot = pA ig ^j^y state invariant under 

rotations around the z-axis, and = ^iCi + ^4C4- this 
choice fulfills the requisite of U{1) symmetry of (g) 
leading to Eq. ([TOj). Referring to the usual Jordan- 
Wigner transformation, we cast Eq. ([1]) in the fermionic 
quadratic form 



(24) 



where {ci,c|} are fermionic operators whose nearest- 
neighbor interaction is described by a (A^+2) x (A^+2) 
tridiagonal mirror-symmetric matrix 



n 



■2/^0 jo 
jo 2h 1 

1 2h 1 



(25) 



1 2h 1 

1 2h jo 
jo 2ho. 



an orthogonal transformation = {Oni} diagonalizes 
H (see Appendix [A]) in terms of Fermi operators 
Cn=^^OniCi and cj which annihilate/create excita- 
tions of energy cj^ [lH, H^- The trivial time-evolution 
of the Cn's entails a time-dependent transformation 



i=0 



where 



(26) 



(27) 



Reminding that ("1,4 = ^(1 ± cr^), ("2 = Cl = ■, a^nd 
that in Eq. © acts on the qubit B, and provided that 



Tii(t) = |w(t)|2+t,(t) 
T44(t) = 1 - V{t) 

T22(t) = -p(<T^+i)«(t)e-W 



where 



U{t) = \UN+lfi{t)\ , 

a{t) = arg[[/Ar+i,o(t)] 



(28) 



(29) 
(30) 



V{t) = \UN+l,N+l{t)f ''"^^' ' " +CN+l{t) , (31) 

AT 

^^w= E u:^m^At)^^^[p'^c]c^'] , (32) 

with p = Tr[Pp^], and P = exp (ztt c|ci) = 

Y[iLi{~^i) is the chain parity operator, which is a con- 
stant of motion. In particular, when p^ is the ground 
state of the chain, p = sign [ det 11] [27^, where O is ob- 
tained from Q by referring only to the chain, i.e., delet- 
ing the first and last row and column. In our case, 
p= (— 1)['^^^°^ ^)/^]^ where [ • ] denotes the integer part, 
meaning that p equals (—1) to the power of the number 
of negative energy fermionic eigenmodes of F. For h = 
it reduces to p= {-l)^^^^l 

Once the above expressions are evaluated, the dynam- 
ics of B directly follows from the initial state of A via 
Eq. (jl]). Similarly, the dynamics of the qubit pair CUB 
follows from the initial state of the qubit pair CUA, via 
Eq. ([7]). Moreover, the time evolution of the magnetiza- 
tion along the chain can be straightforwardly obtained 
from the formalism described above: 

«W) = \U^o{t)\'{a^o) + |^Wi(t)|'(a^+i) + G,{t) , 

(33) 

where G^{t) = 2Q{t) + |/77v+i,o|' + |/7iv+i,iv+i|' - 1. 
The r.h.s. of Eq. (|33|) shows three distinct contribu- 
tions related with the dynamics of the excitations ini- 
tially present in A, B, and the chain F, respectively. For 
h = and N even, Gi(t) = 0, meaning that {(jf{t)) are 
solely determined by two traveling excitations, starting 
from the edges of the chain. These excitations, because 
of the single-particle nature of the Hamiltonian, do not 
scatter, which has great relevance as far as the transport 
properties of the chain are concerned, as discussed in 
Section HVl On the other hand, when the chain is initial- 
ized in the fully polarized state \i) ^^e contribution 
Ci{t) vanishes and Gi{t) has the only effect of redefining 
the range of the magnetization during the dynamics. 



B. Fidelities and concurrence 

We again consider the system ruled by the Hamilto- 
nian ([1]) in the setup described above, so that Eqs. (p!Q|) 
hold. Using Eq. (|28|) we find the entanglement fidelity 
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and the average transmission fidelity of pure states 

J'lfit) = i + lu\t) - ipcosa (a^+i) u{t), (34) 



ent 
Impure / 



T'.Tit) = i + lu^t) - Ipcosa ((T^+i) u{t), (35) 



as well as the concurrence. 

Cc^b"(*) = 

where 



max 



{0,Co} 



(36) 



Co = \p{a%+,)\ u{t) - ^vm-uHt)-v{t)] . (37) 

From the above formulas it appears that the choice of 
the initial state plays an important role [28| : in 
particular, in order to get the largest concurrence it must 
be 



(38) 



meaning that is an eigenstate of P and the qubit B 
is initially in a polarized state, p^ = (i or p^ = (4; as for 
the initial state of the channel, the choices range, for 
example, from its ground state to a fully polarized state. 
Such limitation in the choice of the initial state might be 
overcome by applying a two-qubit encoding and decoding 
on states p^ and p^, respectively [l9, 29]. Similarly, for 
the transmission fidelity the condition 



-p{^N+i) cosa = 1 



must hold. Notice that the l.h.s. of Eq. (|39|) follows from 
the rotation around the z-axis undergone by the state 
during the transmission and can be treated by choosing 
a proper magnetic field ^ or the parity of A^, as well as 
by applying a counter-rotation on the qubit B [28] . 

The above analysis shows that, once condition (|39]) is 
fulfilled the quality of the state and entanglement transfer 
mainly depends on u{t) and increases with it; the residual 
dependence on v{t) suggests that the latter quantity need 
being minimized. In Ref. [28] it has been shown that v{t) 
is exactly zero provided that both V and B are initialized 
in the fully polarized state. On the other hand, since 
Uij{t) entering Eqs. ([29]) and (|3T]) is a unitary matrix, if 
a certain time t* exists such that ii(t*) is very close to 
unity, then v{t*) must be close to zero, no matter the 
initialization. This situation is related with the optimal 
dynamics studied in Ref. [6] : such relation stands on the 
analytical ground developed in the next Section. 



IV. OPTIMAL DYNAMICS 

We now have the tools for determining the conditions 
for a dynamical evolution that corresponds to the best 
quality of the transmission processes. In Appendix [A] 
the algebraic problem of diagonalizing the XX Hamilto- 
nian ([24j) in the case of nonuniform mirror-symmetric 
endpoint interactions, Eq. (|25|) is analytically solved. 
The eigenvalues of Q can be written as 



-h — cos k 



(40) 



in terms of the pseudo-wavevector /c, which takes N-\-2 
discrete values kn in the interval (0,7r): from Eqs. ()A15p 
and (|A16p it follows that these values obey 



7rn + (/Pfc^ 



(n = 1, 7V+2) , (41) 



with 



/ oot k \ 

ipk = 2k-2cot-^i^-^j G(-7r,7r), (42) 



■fo 



(43) 



where we have set ho = h. From the above equations 
it follows that the /c's correspond to the equispaced val- 
ues nn/ (A^+3), slightly shifted towards 7r/2 of a quantity 
which is smaller than tt/ (A/'+3), so that their order is pre- 
served: therefore k can be used as an alternative index 
for n, understanding that it takes the values /c^, as done 
in Eq. (|4Q]) . According to the conclusions of the previous 
section, we focus on the transition amplitude, Eqs. (|27|) 
and (f29|) - (f3Q|) . which explicitly reads: 



(39) UN+iMt) = w(t)e^"^*^ = - E '^('^") e 



z(7rn— ojfc^ t) 



(44) 



where, after Eq. ()A2ip . it is 



p{k) 



1 



^_ A(l+A) 
A2 + cot2A; 



(45) 



and mirror symmetry is exploited according to Eq. ()A3p : 
the transition amplitude above is a superposition of phase 
factors with normalized weights, p{kn) = 1, entailing 
\u{t) \ < 1 with equality holding when all phases are equal. 
The distribution p{k) is peaked at k = ko = 7r/2 and its 
width is characterized by the parameter A, Eq. (|43|) . so 
that the smaller jo the narrower p{k). 

As u{t) essentially measures the state-transfer qual- 
ity, the condition for maximizing it at some time t*, i.e. 
u{t^) 1, is that all phases irn — ook^f" almost equal each 
other. Assume for a moment that the /c's be equispaced 
values, as in (jA8p . and that the dispersion relation be 
linear, uok = vk: then Eq. (jHj) would read 



u(t) 



^Z7rn{l — t/t*) 



(46) 



with r=(A/'+3)/v, so that ^x(r) = En/^(^n) = l, i.e., ah 
modes give a coherent contribution and entail perfect 
transfer. On the other hand, in our case is nonlin- 
ear in /c, and the kn are not equally spaced due to the 
phase shifts (|42]) entering Eq. (|4T]) . so generally the dif- 
ferent modes undergo dispersion and lose coherence. 
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A. Transfer regimes 

The dependence of A upon jo reveals the possibihty of 
identifying different dynamical regimes, characterized by 
a qualitatively different distribution Eq. (|45|) . and 

hence, as for the transfer processes, a different behavior 
of the transition amplitude u{t). For extremely small jo 
the distribution p{k) can be so thin that (for even N) 
only two opposite small eigenvalues come into play, say 
differing by 5uj^ and perfect transmission will be attained 
at a large time t = 7r{Suj)~^ (for odd TV there is a third 
vanishing eigenvalue at /c = 7r/2 and still the two identical 
spacings Su do matter). This is the Rabi-like regime also 
mentioned in the Introduction. 

A different regime is observed when jo is increased: 
a few more eigenvalues come into play and it may oc- 
cur, in a seemingly random way, that their spacings be 
(almost) commensurate with each other, i.e., they can be 
approximated as fractions with the same denominator 
yielding phase coherence at = ttK. By recording the 
maximum of u{t) over a fixed large time interval T, as jo 
is varied (see Ref. [sl]), a rapid and chaotic variation is 
observed. This regime is clearly useless for the purpose 
of quantum communication. 

As jo further increases, the ballistic regime eventually 
manifests itself: p{k) involves so many modes that com- 
mensur ability is practically impossible, and a more regu- 
lar behavior with short transmission time ^ N sets in. 
The ballistic regime is characterized by relatively large 
values of ii(t*. A) which is the quantity plotted in Fig. [21 
reporting numerical results for increasing chain lengths. 
It appears that each curve shows a maximum for a par- 
ticular optimal value of A = A°p*(A/') or, equivalently, of 
jo = jQ^*(7V): such maxima are remarkably stable for very 
high N and yield very high transmission quality. In Ta- 
ble [J we report some of the optimal values A°p* (N) and 
j°P*(7V) for a wide interval of chain lengths. This last 
'ballistic-transfer' regime is the one we are interested in, 
since it has two strong advantages: first, the transmis- 
sion time ^ N is the shortest attainable, and second, 
the maximum value i^(t*,A°P*) of 'u(t*,A) is such that 
one can achieve very good state transfer, e.g., the corre- 
sponding transmission fidelity is far beyond the classical 
threshold, even for very long chains. 

The above analysis gives a physical interpretation of 
what is observed in Fig. 3 of Ref. [3], where the Rabi- 
like, intermediate and ballistic regimes emerge. 

A qualitative picture of the ballistic regime can be 
obtained by viewing the transition amplitude (|44|) as 
a wavepacket with components. It can be evalu- 

ated by progressively adding the contributions from sym- 
metric eigenvalues, i.e., for odd N summing between 
(Ar+1)/2T^, for ^ = 0, 1, (7V+l)/2. This yields the 
partial sum ui(t*) shown in Fig. [3l together with the 
corresponding frequency and density. One can see that 
the amplitude increases only over the modes of the linear- 
frequency zone, i.e. where frequencies are equally spaced, 
indicating that only those wavepacket components whose 
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FIG. 2. (color online) Value of i^(t*, A) as a function of A, 
for different wire lengths N . t* is obtained numerically by 
maximizing Eq. (|44|) around t:^ 




FIG. 3. (color online) Partial sum of the amplitude Ui{t*) vs i 
for + 2 = 51 and jo = 0.58, together with the corresponding 
frequency and density. 



frequency lies in such zone play a role in the transmission 
process. 



B. Ballistic regime and optimal values 

From the above reasoning, since the modes contribut- 
ing to the amplitude lie in a range of size A around /cq, 
in order to get high-quality transfer processes it is neces- 
sary that the corresponding frequencies be almost equally 
spaced, meaning that is approximately linear in n. 
Actually, ujk has an inflection point in ko: its nonlinearity 
is of the third order in k — ko and the modes close to ko 
satisfy the required condition. However, from the phase- 
shifts (|42]) a further cubic term arises, which depends on 
A. As A varies with jo, the latter can be chosen so as to 
eliminate the cubic terms, yielding a wide interval with 
almost constant frequency spacing. The latter can be 
expressed just as the derivative of with respect to 
n, dnOJkn = sin/c dnk. The last term is evaluated from 
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N + 2 


25 


51 


101 


251 


501 


1001 


2501 


5001 


10001 


25001 


50001 


100001 


250001 


500001 


^opt 


0.245 


0.182 


0.139 


0.098 


0.075 


0.058 


0.042 


0.033 


0.026 


0.019 


0.015 


0.012 


0.009 


0.007 


Jo 


0.628 


0.556 


0.494 


0.422 


0.374 


0.332 


0.284 


0.252 


0.224 


0.192 


0.171 


0.152 


0.130 


0.116 




0.972 


0.953 


0.936 


0.916 


0.902 


0.891 


0.880 


0.873 


0.868 


0.862 


0.859 


0.857 


0.854 


0.853 



TABLE I. Optimal values A°p' and the corresponding and w(t*, A°p') (see text for details), for different N. 



Eqs. (j4T|) and (|42l). 



1-A 



2(l-A^)cos^fc 
A[A2 + (1-A2) cos2 k] 



(47) 
(48) 



SO that 



sin k 



TT 



1-A2 1 



where T = A/'+3 + 2 (1-A)/A is the arrival time. It fol- 
lows that one can minimize the nonlinearity of uJkr^ by 
setting the width to the value Aq satisfying 



Ao 



(1-A^) 



.1/3 



22/3JY-1/3 



(50) 



and jo — 2^/^A^~^/^ for large N . Therefore the main 
mechanism that produces an optimal ballistic transmis- 
sion is that of varying the endpoint exchange parameter 
to the value jo that 'linearizes' the dispersion relation. 
Actually, if the corresponding Ao = A(jo) is such that 
p(k) exceeds the region of hnearity, further gain arises by 
lowering jo so as to tighten the relevant modes towards 
/co- However, at the same time cj^^ becomes less linear 
and the trade-off between these two effects explains why 
a maximum is observed. This is well apparent in FiglH 
where for different values of A the shapes of dnOOk can be 
compared with the excitation density p{k): for A = Ao 
the density still has important wings in the nonlinear 
zone, so the optimal value A°p* turns out to be smaller. 

The dynamics in the ballistic regime is best illustrated 
by the time evolution of the magnetization Eq. (|33]) along 
the chain, plotted in Fig. [5] when the initial state is 
It) ^ lii • • • i) ^ It)- The initial magnetizations at the 
endpoint s generate two travehng wavepackets: for non- 
optimal couplings (jo = 1, upper panel) they change their 
shape and quickly straggle along the chain; for optimal 
coupling (jo = jo , lower panel) they travel with minimal 
dispersion. This confirms that the coherence is best pre- 
served when the optimal ballistic dynamics is induced: 
In the next session we show that to such dynamics do in 
fact correspond high values of the quality estimators for 
the state and entanglement transfer. 



1.0 



0.9 



0.8 



0.7 



1 + (2 7^-2) co^'^ + O(cos^fc)] , (49) 7T/4 





N+2=51 


V ->^'^ \ X 

k \ A"?'/^ 




/ '0.1/ 






p(k) ^ 
p(TT/2) : 



1.0 



0.5 



71/2 



3Tr/4 



0.0 



FIG. 4. (color online) 'Group velocity' Vk = [{N-\-3)/7t] dnOOk^ 
and p{k) vs k for different values of A. The thicker curves 
correspond to Ao =0.3944 (|50|) that gives the flat behavior at 
/Co, and to A°P':^ 0.1825. 



V. INFORMATION TRANSMISSION 
EXPLOITING OPTIMAL DYNAMICS 

The requirement Eq. (|39|) . means that the state is 
not rotated by the dynamics when it arrives on site B, 
though during the evolution it may undergo a rotation 
around the z axis. In Ref. [20] it has been shown that 
a= — ^{N^l) at the transmission time t*. Therefore, 
also without applying a counter-rotation on qubit B [28| , 
condition (|39|) can be fulfilled by choosing A^ = 4M±1 
where the sign ± is given by (|38|) and thus depends on 
the initial state of the chain. In the following we assume 
that conditions (|38|) and (|39|) are always satisfied. 

Let us consider for the moment that F and B are ini- 
tially in the fully polarized state |tt • • • t) It)- In that 
case v{t) = {) and the transmission fidelities and (|35]) , 
as well as the concurrence (|36|) , only depend on, and 
monotonically increase with, u{t). The best attainable 
information transfer quality corresponds therefore to the 
maximum amplitude iiopt = 'w(t*, A°p*). In Fig. [6] and in 
Table |I] we report these values together with the corre- 
sponding optimal A°P* as a function of the chain length 
A/" in a logarithmic scale; the inset shows that A°p* obeys 
the same power-law behavior predicted in Eq. (|5Q]) for 
Ao. Fig. [6] also shows that for larger and larger N the 
maximal amplitude iiopt does not decrease towards zero, 
but it rather tends to a constant value of about 0.85, 
which is surprisingly high, as, e.g., it corresponds to an 
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1.00 







100 200 300 400 , 500 



FIG. 5. (color online) Dynamics of the magnetization af{t) 
at time t and site i when a) jo — 1 and b) jo — 3q^^ • The initial 
state of the whole system is |t) |ii • • • |t) "the length 
of the chain is + 2 = 250. 
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FIG. 6. (color online) Behavior of the maximum attainable 
amplitude i^opt and (inset) of the corresponding optimal value 
of A°^^ vs logarithm of the chain length N . The horizontal 
dashed line is the infinite N limit of Uo^t- 
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average fidelity ^ab(^*)^0-9- This is indeed true: We 
show in Appendix [B] that in the limit oi N ^oo the op- 
timized amplitude tends to i^opt = 0.8469. Basically, this 
tells us that it is possible to transmit quantum states 
with very good quality also over macroscopic distances. 
From Eq. (|BT4|) we can derive the asymptotic behavior 
of the optimal coupHng 

- 1.030 iV-^/^ . (51) 

In the optimal ballistic case the channel initialization is 
not crucial, as different initial states satisfying (|38|) give 
rise to almost the same dynamics as discussed at the end 
of subsection III B. In fact, the term CN+iit) entering 
Eq. ([3T]) essentially embodies the effect of channel initial- 
ization and it is expected to be small at t = t*. This is 
apparent in Fig. [71 where for jo = Jo^*, CN+i{t*) stays 
well below 0.1 for TV as long as 1000. 

The transmitted entanglement, as measured by the 
concurrence (|36|) . is shown in Fig. [8] as a function of jo 
and t, with the channel initially prepared in its ground 
state. As expected, the peak of the transmitted con- 
currence is observed for jo = jo^*; away from jg^* the 
quality of transmission falls down because u{t*) decreases 



FIG. 7. (color online) CAr+i(t*) for different initial states of 
the chain (ground state, ant i- ferromagnetic Neel state, and 
series of singlets [28|] ) when jo = jo and jo = 1 • The results 
for a series of singlets are numerically indistinguishable from 
those with the Neel state. 




500 



0-4 0.6 0.8 1 

jo 



FIG. 8. (color online) Evolution of the concurrence Oq-q vs 
jO and t. The length of the chain is + 2 = 250. 
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FIG. 9. (color online) Minimum fidelity vs time for differ- 
ent values of jo. The length of the chain is N -\-2 = 251 and 
0.422. 

and, accordingly, v{t*) is allowed to increase. In fact, in 
the non-optimal ballistic case the quality of entanglement 
transfer does depend on the the initial state of the chan- 
nel [m, ^0]; for instance, when jo = 1 and the chain is 
initially in its ground state, the contribution of the over- 
lap terms Tr[p^cjcj/] in Eq. (|32]) is not quenched by the 
dynamical prefactors, and higher values of CAr+i(t*) (see 
Fig. [7j) inhibit the transmission of entanglement even if 

7i(r) ^ 0. 

The effect of the optimization of jo is clearly evident 
in the time behavior of the minimum fidelity, Eq. (j22j), 
reported in Fig. [9l The peak of J^^^{t) for Jq^* occurs 
at the arrival time + s with a time delay s that 

agrees with the asymptotic value s 2.29 A/"^/^ derived in 
Appendix [Bl The 'reading time', i.e., the time interval 
during which the qubit B keeps being in the transferred 
quantum state, is A~^, as the same figure also 

shows; note that, in the optimal case, increases with 
N according to the asymptotic behavior ^r 1.89 N^^^. 



VI. CONCLUSIONS 

In this paper we have shown that high-quality quan- 
tum state and entanglement transfer between two qubits 
A and B is obtained through a uniform XX channel of 
arbitrary length by a proper choice of the interaction 
jo between the channel and the qubits. The value of such 
interaction is found to control the transfer regime of the 
channel, which varies, as jo increases, from the Rabi-like 
one, characterized by very long transmission time, to an 
intermediate regime, which turns useless for the purpose 
of quantum communication, and finally becomes ballistic 
for jo of the order of the intrachannel interaction. 

In order to get coherent transfer in the ballistic 
regime, it is desirable that the /c-density of the travel- 
ing wavepacket generated by Alice's initialized qubit A 
be narrow and concentrated in the linear zone of the dis- 



persion relation, i.e., with equispaced frequencies. As the 
parameter jo controls both the width of the /c-density and 
the spacings of the frequencies entering the dynamics, one 
can therefore improve the transmission quality up to a 
best trade-off arising for an optimal value jopt (^) which 
for large N behaves as jopt (7V)c^l.03 7V-i/^. Remark- 
ably, we have found that for such a choice the quantum- 
state-transfer quality indicators are very high and, in- 
deed, have a lower bound for ^ oo that still allows to 
efficiently perform quantum-information tasks: e.g., the 
average fidelity of state transmission is larger than 90 %. 

The ballistic regime ensures fast transmission on a time 
scale of the order of A^, at variance with the Rabi-like 
regime, and in the optimal case the reading time increases 
as A^^/^. It is also to be noted that, if experimental set- 
tings constrain to a given value Jq^^, yet one can opti- 
mize the chain length in such a way that Jq^^ = jQ^^{N). 
The only requirement on the initial state of the receiving 
qubit B and the spin bus is to possess U{1) symmetry, a 
condition that can be fulfilled by several configurations 
concerning the spin bus, ranging from the fully polarized 
state to the highly-entangled ground state. If a large 
magnetic field can be switched on during the initializa- 
tion procedure (in order to fully polarize the channel), 
and switched off as soon as the transmission starts, then, 
from our analytical treatment it emerges that tempera- 
ture is not a major issue as far as the dynamical evolution 
of the channel is concerned, though low temperatures are 
obviously necessary to protect the qubits from phase and 
amplitude damping due to the solid-state environment. 
To judge if the proposed scheme identifies a reasonable 
experimental framework, let us estimate the magnitude 
of the involved physical quantities. Consider a solid- 
state implementation with lattice spacing of about 10 A 
and intrachain exchange J :^ 10^ K. A quantum state will 
then be transferred with fidelity 90 % along a channel of 
length 1 cm (A^::^10^) using jo ^1.03 N'^^^ J c^7.0 K, 
with transmission time t = N h/ {kBJ) —0.7b fis and read- 
ing time tR0^1.9N^^^ ^/(/cb J) ^ 0.03 ns. 
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Appendix A: Quasi-uniform tridiagonal matrices 

The matrix Eq. ([25]) . can be written diS Q = —h 
\M where 



X y 
y 1 
1 1 



(Al) 



1 1 
1 y 

y X 



is a square tridiagonal matrix of dimension M = A/" + 2, 
and X = 2{ho — h) and y = jo- This real symmetric 
matrix is diagonalized by an orthogonal matrix 0{x^y)^ 



M 



Afe ^kk' , 



(A2) 



and it is known that {%) if ?/^0 the eigenvalues are 
nondegenerate [31], (n) the eigenvectors corresponding 
to the eigenvalues ordered in descending order are alter- 
nately symmetric and skew symmetric [32], i.e.. 



(A3) 



The eigenvalues are the roots of the associated charac- 
teristic polynomial 



Xm(A;x,?/) =det[\- M{x,y)\ . 



(A4) 



In the fully uniform case the characteristic polynomial 
is r]^{X) = x„(A;0, 1) and one easily obtains the recur- 
sion relation 



(A5) 



that can be solved in terms of Chebyshev polynomials of 
the second kind. 



where 



sin(M+l)^ 
sin k 



A = 2 cos k , 



(A6) 



(A7) 



SO the eigenvalues of A^(0,1) correspond to M discrete 
values of /c. 



TT n 



M+1 

the corresponding eigenvectors are 



[n = l,...,M] 



Ofcz(0,l) = 



M+l 



sinki . 



(A8) 



(A9) 



The general determinant ()A4p can be expressed in 
terms of the t^^'s by expanding it in the first and then 
in the last column, 

Xm = {X^-2xXW) r]^_^-2y^{X-x) r]^_^ + / r]^_^ , 

(AlO) 

and using Eq. (|A5p one can eliminate the explicit appear- 
ances of A, 

Xm ~ Vm ~ Vm-1 ~^ X Vm-2 

+(lV)[2r?M_,-2x7?^_3+(lV)^M-4](All) 

By rewriting Eq. (|A6p as sink r]^ = Im ^e**^^+-^^^] and 
defining 

z^ = l-y\ zl = z^e-^'^, Xk=xe-'^ , (A12) 



Uk = I — Xk ^ z1 = e { [{2—y'^) cos k — x] ^ i y^ smk^ , 

(A13) 



Eq. (lAlip takes the form 

sink x^(i^) =Im{e^(^+i)^ ul] . (A14) 

The secular equation Im{e*^^+^^^ = entails that 
when k corresponds to an eigenvalue the quantity in 
braces is real and equal to either ±.\uk\^] by Eq. (|A13p it 
turns into sin [(M+l)/c — Lpk\ = 0, with the phase shift 



ifk = 2k — 2 tan 



y'^ sin k 



(2—y'^) cos k — X 
so the M eigenvalues correspond to 

7rn + ipkr^ 



kn 



M+l 



(n = l,...,M) . 



(A15) 



(A16) 



We are interested in the squared components of the 
first column of the diagonalizing matrix which can 
be expressed as [31] 

Ol^ = = -'"f , (A17) 



^AXm(A/c) 



dkXMik) 



where k assumes the values ()A16p and (A; x, ?/) is the 
characteristic polynomial associated to the first minor 
matrix A^n that, expanded in the last column and using 
Eq. (|A5]) . reads 

=det[A-Mii(x,^)] 



= {^-X)Vm-2 -y Vm-3 



MS . (A18) 
(A19) 



Then the numerator of Eq. (|A17p is 

sink =lm{e'^^Uk} , 

while from Eq. ()A14p one has 

sinfc dkXMik) = (M+l) Re{eWi)fe ^1} 

+2 Imje^^^+i^'^ Ufcw'J , (A20) 
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where the argument of Re is real indeed; retaining only 
the dominant term for M ^ 1 Eq. (|A17p becomes 



M+1 [(2-^2)cosA:-x]2+^4gin2^ 



. (A21) 



For X = the above expression is in agreement with 
Ref. [3]. In the most common case x<2—y'^ the maxi- 
mum /co of Of^^ is located at 



cos /co 



(A22) 



so that switching x on the maximum shifts from 7r/2 to 

fco = I - sin-1 ^ ; (A23) 

the 'eigenvalue' corresponding to the maximum is 

2x 



Ao = 2 cos ko = 



(A24) 



so that for ~ 1, the maximum shifts the 'energy' lin- 
early with X. Expanding around the maximum, the 
leading behavior is found to be a Lorentzian, 



Oi 



y 



~ M+l + [(2-^2)2 _ ^2](^ _ ^q)2 

whose width (HWHM) is given by 



, (A25) 



(A26) 



When X and y are small, ko {tt—x)/2 and A ~y^/2, so 
X rules the position of the peak, while y determines its 
width. 



Appendix B: large- limit of the amplitude 

The transition amplitude, Eq. (ji^ . in the case of odd 
Ar = 2M-l reads 



u{t) 



A(l + A) 



^ ^t{7nn — t sin Qm) 



E 



+ 3 A2 + tan2 ' 

171 = — M 



(Bl) 



where the summation has been made symmetric through 
the change of variable q = 7r/2 — k. The shift equa- 
tion (|4T]) turns into 



TT m = (7V+3) Qm + TTCPq^ 



with 



TTcpq = 2^ tan" 



tang 



(B2) 



(B3) 



In the limit N ^ 00 one can write the sum as an integral 
setting 



N^3 



TT V N+3J 



As we deal within the region of the optimal value of 
A A/"-^/^ 0, we have 



{t) = lim A 



2[(A/'+3)g+7r(^q— t sin q] 

A2 + tan^ 5 



(B5) 



Writing the arrival time as t = N-\-3-\-s^ where s is the 
arrival delay, one has then 



lim A 

t^oo 



f fj^q ^t[t{q- sin q)-sq-\-7ripq] 



n dq 



A2 + tan^ q 



(B6) 



The relevant q^s are of the order of A^N~'^^^ 0, so 
we change to g = Ax, with x of the order of 1; keeping 
the leading terms for A ^ 0, 

t {q - sin q) — > x^ , (B7) 



TT 

A dq 
A2 + tan^ q 



6 

2 tan~^ X , 
dx 



(B8) 
(B9) 



and defining the rescaled counterparts of the arrival time 
t :^ A/" and of the delay s - A^^/^, 



A^ 

6 



A 5 ; 



the final asymptotic expression results 



(BIO) 



(Bll) 



that can also be rewritten in the form of a simple summa- 
tion of phases by introducing the variable z= tan~^ x. 



Jo 



cos(rtan^ z — (T tan 2^ + 22;) . (B12) 



As in the finite- A^ case, one has to maximize Uoo{t, a) by 
finding the optimal values of a and r. For r = it is easy 
to evaluate Eq. ()Blip analytically. 



^00(0, cr) =2e~ 



cr ; 



(B13) 



it is maximal for cr = 1, giving u{0, 1) = 2e~^ ^ 0.736, to 
be regarded as a lower bound to the overall maximum 
of Uoo{T^a). The overall maximization has been per- 
formed numerically using Eq. ()B12p . It turns out that 
the maximum corresponds to a = 1.2152 and r = 0.02483, 
and amounts to iioo(0.02483, 1.2152) =0.84690, in agree- 
ment with the behavior shown in Fig. [6l The resulting 
scaling, from Eq. (jBlOp . tells that asymptotically 



A 0.530 A"-^/^ , s 2.29 N^^^ . 



(B14) 
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